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1 Introduction 



Heavy flavour decays play an outstanding role in the determination of fundamental parame- 
ters of the standard model and for the development of a deeper understanding of the QCD 
dynamics. In order to achieve the ambitious scientific goals it is important to combine the 
experimental and theoretical efforts. However, this is not always easy. For example, while 
the theory of inclusive decays is most advanced, inclusive measurements are generally quite 
difficult. Conversely, exclusive decays into few-body final states are often much easier to detect 
and reconstruct, but the theoretical methods to treat exclusive processes are still rudimentary 
when hadrons are involved. In view of the exciting experimental prospects at future bottom 
and charm factories, where many new exclusive channels will become accessible, it appears 
particularly desirable to have a reliable quantitative theory. 

As well known, the main problem is the complicated interplay of weak and strong interac- 
tions, and most of all the influence of long-distance dynamics. Clearly, a complete theoretical 
understanding of the decays of hadrons containing a heavy quark is only possible in a nonper- 
turbative framework. This is most obvious for exclusive nonleptonic decays and for processes 
involving strong couplings among hadrons such as B*Btt or D*Dtc. Yet, what appears as 
a problem for determining electroweak parameters from exclusive measurements can also be 
regarded as a chance to improve and test our understanding of hadronic physics in QCD. 

The nonperturbative approaches applied in this field include lattice approximations, heavy 
quark effective theory, chiral perturbation theory, l/iV c -expansion, QCD sum rules, and quark 
models. Each of these methods has advantages and disadvantages. For example, quark models 
are easy to use and very good for intuition, however, their relation to QCD is unclear. In the 
other extreme, lattice calculations are rigorous from the point of view of QCD, but they can (at 
least presently) be applied only to relatively simple matrix elements related to decay constants 
and form factors. Moreover, there are uncertainties connected with the necessary extrapola- 
tions to the physical masses of the light (u,d) and heavy (b) quarks. Also the application of 
effective theories is restricted to a rather narrow class of problems, and the corrections are 
often substantial. More importantly, in many cases the latter cannot be calculated directly 
within the same theory, and therefore bring in model-dependence after all. 

QCD sum rules |l[ represent no exception is this respect. On the one hand, they are based 
on first principles. On the other hand, one has to introduce new elements such as vacuum 
condensates or light-cone wave functions which cannot (yet) be calculated directly in QCD. 
However, since this nonperturbative input is of universal nature, sum rules nevertheless possess 
a high predictive power. Furthermore, the sum rule approach is particularly well suited for 
heavy quark physics, since the approximations which have to be made in practice can be well 
justified when a heavy mass scale is present. Last but not least, the sum rule approach is 
less limited in applicability than some of the other methods mentioned, or the limitations 
are different. In this paper we will exemplify some recent applications of the QCD sum rule 
technique to B and, for comparison, also to D decays. 

The simplest and therefore most instructive example for introduction is actually provided 
by the calculation of heavy meson decay constants. They are defined by matrix elements of 
quark currents, for example 

m b (B | 6ry 5 d | 0) = f B m 2 B , (1) 

and determine leptonic decays such as B — ► tv t . In addition they enter in the factorizable 
amplitudes of nonleptonic decays. Without going into details, which for example can be found 
in 0, we only mention that the QCD sum rule approach to (|I]) is based on an analysis of 
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the two-point correlation function (0 | T{b(x)ij 5 d(x), J(0)i7 5 6(0)} | 0) using operator product 
expansion, perturbation theory, and dispersion relations. T in the above current-correlator 
denotes the time-ordering operation. 

For brevity, we directly proceed to more complicated cases and begin with the calculation 
of form factors. In section 2, we describe how sum rules on the light-cone can be used to 
derive the B — > 7r and D it form factors. Comparison with recent experimental results on 
CKM-suppressed B and D decays as well as the determination of \V ub \ are discussed in section 
3. Finally, section 4 is devoted to the main theoretical problem in heavy quark decays that 
is the complete calculation of exclusive nonleptonic amplitudes. We briefly demonstrate the 
failure of strict factorization, the usual phenomenological procedure, and a first attempt to 
estimate the nonfactorizable piece of the amplitude for B — > J/ipK. 



2 The B — » 7T and D — > 7r form factors 

Recently, the CLEO collaboration [^] announced the first measurement of the decay B — > tiIu. 
This process plays a very important role for the determination of the CKM parameter V u b. 
The decay amplitude is completely determined by the form factor f^ip 2 ) entering the matrix 
element 

<7r(q)\u^b\B(p + q) >=2f+(p 2 )q,+ [fUp 2 )+f B (p 2 )}p„ (2) 

where the assignment of the momenta p and q is quite obvious. In |3J] this form factor was 
calculated from a QCD sum rule on the light-cone. In the following we briefly outline the 
method. 

One starts from the correlation function of the relevant vector and pseudoscalar currents 
between the vacuum and the on-shell pion state: 



*V,(P + <!?)= • (4) 

. I ty-. 1^ J-* \A \ 



Fp(p,q)=i J d A xe ipx (v{q) | T{u(x) lt ,b(x), b(0)i l5 d(0)} | 0) 

= F(p 2 , (p + qf)q^ + F(p 2 , (p + q) 2 ) Pfl . (3) 

Writing a dispersion relation in the variable (p + q) 2 , one readily finds the B meson ground 
state contribution to the invariant function F(p 2 , (p + q) 2 ) by inserting a complete sum of 
states between the currents in ((^), focusing on the term \B)(B\ and using the relations (fj) 
and (g): 

2m 2 B f B fUp 2 ' 
m b (m 2 B - (p + g) 5 

The ellipses in the above denote the contributions from the excited B and from continuum 
states. 

The main idea of the method is to apply the operator product expansion (OPE) to the 
product of currents in the correlation function (§) at spacelike momentum (p + q) 2 < such 
that the b-quark propagating from point x to point is highly virtual. Moreover, the expansion 
is performed near the light-cone, x 2 = , rather than around x = 0. Consequently, one ends 
up with a series of bilocal rather than local operators. More definitely, after contracting the b- 
quark fields in (|3]) and expanding the 6-quark propagator near x 2 = 0, the correlation function 
can be expressed in terms of quark-antiquark and quark-antiquark-gluon matrix elements of 
the type 

(ir(q) | u(x)T a d(0) | 0) and (vr(g) | u(x)g s G ,J ' u (vx)r b d(0) | 0) , (5) 
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respectively, where T a ^ denote certain combinations of Dirac matrices and the auxiliary vari- 
able v varies from to 1. These matrix elements are essentially nonperturbative objects. 
Most importantly, they are universal and process- independent. In particular, the same matrix 
elements enter in the corresponding expansion of the correlation function interpolating the 
D — > 7r form factor. The latter can be immediately obtained from (|3|) by replacing b by c and 
B by D. 

The vacuum-to-pion matrix elements introduced above can be parameterized in terms of 
so-called light-cone wave functions with given twist || |6], [7]. The most important of them is 
the leading twist 2 wave function (p^u, f/,) defined by 

(7r(g)|u(x)7 M 7 5 d(0)|0) = -iq^U I due iuqx ^ n (u, fj) + ... , (6) 

J o 

where /z -1 is a characteristic distance scale and u is the fraction of the pion light-cone mo- 
mentum q + g 3 carried by the constituent quark. The ellipses in (||) stand for the higher 
twist components. Asymptotically, that is at \i — * oo, QCD perturbation theory implies 
<p n (u, oo) = 6u(l —u). However, at the physical scale \i ~ m&, at which the OPE is applied to 
the correlation function (|3]), important nonasymptotic effects are to be expected. These are 
parameterized by the coefficients aj(/i) in the following series of Gegenbauer polynomials: 

<p v {u,fi) = Gu(l-u)[l+ £ a i (^)Cf /2 (2u-l)] . (7) 

i=2,4,.. 

Although the input values of the coefficients at a fixed low momentum scale \x are 

unknown, the scale-dependence of eij(//) is dictated by the renormalization group. From what 
is said above, it is also clear that Oi(/x) —>■ for \i — > oo. Over the years a great deal has been 
learned about these wave functions. They have been classified up to twist 4 and the asymptotic 
form has been determined. Also various nonasymptotic corrections have been estimated from 
QCD sum rules for light hadrons. In our calculation we have used the wave functions given 
in 0. 

Substitution of (^) into the light-cone expansion of the correlation function (^) yields a 
representation of the invariant amplitude F in terms of perturbative and nonperturbative 
parameters and functions of QCD. The leading term of it reads 

F(p\(p + g ?)=m b ,J 1 f +■■■■ (8) 

Jo mi — [p + uqy 

Equating now the QCD result (§) in the region of validity at (p + q) 2 < with the dispersion 
relation (^) one obtains a raw sum rule relation for /b/b(p 2 ). The rest of the calculation 
follows the usual QCD sum rule procedure: Borel transformation in (p + q) 2 and subtraction 
of the contribution from higher states invoking quark-hadron duality. Details can be found 
in the original paper 0. One finally arrives at an expression for the desired form factor, the 
dominant term of which is given by 

f+r 2x f^ m b [ ldu ( ( \ i \ ( m B m 2 b -p 2 (l-u)\ 

Here, M is the Borel mass parameter, and the scale /i& is of order of the characteristic virtuality 
of the correlation function, [i 2 = m 2 B — m 2 . The integration limit A = (m 2 — p 2 )/{sq — p 2 ) 
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depends on the effective threshold So above which the contribution from higher states to the 
dispersion relation (^) are cancelled against the corresponding piece in the QCD representation 
(H). The twist 3 and 4 contributions denoted by ellipses have also been calculated in [|], ||, 
but are not shown here for brevity. They involve further quark-antiquark wave functions as 
well as three-particle quark- ant iquark-gluon wave functions. 

An important detail which should be mentioned here concerns the treatment of f B in 
the sum rule (|9]). For consistency, f B is replaced by the appropriate two-point sum rule, 
disregarding the 0(a s ) gluon corrections since the latter are also not yet included in (H). 
This procedure considerably decreases the sensitivity of f B to the choice of rrib and sq as 
demonstrated in H, |§. Furthermore, we observe a remarkable stability of (|9]) with respect to 
a variation of the Borel parameter M. The appropriate range of M 2 is found by requiring 
the higher twist terms and simultaneously, the contributions from heavier states to remain 
subleading. 

In Fig. la the form factor f B {p 2 ) is plotted versus p 2 for the central value M 2 = 10 GeV 2 
of the fiducial range. On general grounds, one can expect the sum rule (|9|) to be valid up to 
a timelike momentum transfer p 2 = ml — 0(1 GeV 2 ). Actually, the calculation shows that 
above p 2 ~ 17 -j- 20 GeV 2 the stability in M 2 is lost. In Fig. 2 we compare our prediction 
with the results obtained in a quark model || and from a three-point sum rule involving local 
quark and gluon condensates instead of light-cone wave functions |IJ . 



Obvious replacements in (Q) and ([5]) yield the corresponding sum rule for the D — > n 
form factor. The numerical result for /-(p 2 ) is plotted in Fig. lb taking M 2 = 4 GeV 2 in 
accordance with the reliability criteria pointed out above. Note again that this calculation 
should only be trusted at p 2 < m 2 — 0(lGeV 2 ). 

For easy comparison with other theoretical results (a rather complete compilation can be 
found in JTTJ] ) we also quote our predictions at p 2 = 0: 



/^(0) = 0.29 ±0.01 , fS(0) = 0.66 ±0.03 . (10) 

The uncertainty above only refers to the variation with M 2 within the fiducial range. We have 
also investigated the sensitivity to nonasymptotic effects in the light-cone wave functions. To 
this end we have recalculated f B and f£ using purely asymptotic wave functions and find that 
the form factors change by less than 10 % . This gives additional confidence in the results 
exhibited in Fig. 1. 

For an accurate theoretical analysis of the semileptonic decays B — > nlv and D — > nlu one 
also has to know the form factors at larger values of p 2 , ideally up to the zero recoil points 
[m B — m w ) 2 and {m D — m„-) 2 , respectively. At such large p 2 the form factors are expected to 
be dominated by the B* and D* poles, respectively. In this approximation, 

fW) = 9 f ?: 9B '?, 2 , , (ii) 

where Qb'Btt is the B*Bn coupling defined by 

(B*-(p)n + (q) | B°(p + q)) = -g B * B7T q^ . (12) 

It is one of the advantages of the light-cone sum rule method that the B*Bn coupling can 
be calculated || from the same correlation function @. One only has to set up sum rules with 
the respect to both the B- and B*- channel considering a dispersion relation in (p + q) 2 and 
p 2 , simultaneously. After double Borel transformation of the invariant amplitude (§) involving 
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the mass parameters Mi and M 2 one obtains an expression for the product f b$ b*9b*b-k the 
leading twist term of which depends on the pion wave function <p n (u, fi) at the fixed momentum 
fraction u = M 1 2 /(M 1 2 + Mf ) ~ 1/2: 



f sf B*9B* Bit 



m 2 h U m l +m B* 
— e 2Af 2 



m B mB* 



M 2 [e~^ - e"3?r]^(l/2) + . . . 



(13) 



with M 2 = M 2 M 2 /(M 2 + Mf). Taking <p v (l/2,fi = 0.5 GeV) = 1.2 ± 0.2 together with 
the corresponding values of the higher twist wave functions and dividing by the values of 
and fs* as given by the appropriate two-point sum rules, we find, numerically, 



9B* Bit = 29 ± 3 



(14) 



where the indicated error again quantifies the variation within the fiducial range of the Borel 
parameters. 




2.5 3 

p 2 [GeV 2 } 

Fig. 1. The form factors for the transitions (a) B — > n and (b) D — > it as predicted by light-cone 
sum rules (solid) in comparison to the single-pole approximation (dashed) with the normalization 
constants gB*B-w an d 9d*d-k, respectively, calculated by the same method. 



Expressions analogous to flllD and (12 ) hold for fp(p 2 ) in the region near p 2 ~ (m^ — m^) 2 . 
Following the same procedure as the one described above, we get 



g D * m = 12.5 ± 1.0 . 



(15) 
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Fig. 2. Comparison of the B — > tt form factor as derived from the light-cone sum rule (solid) with 
the predictions of a three-point sum rule jYZ|/ (dashed) and a quark model j^j (dash- dotted). 



which is consistent with the limit derived from recent ACCMOR \I2\ and CLEO |TB|] mea- 
surements. 



The single-pole approximation for the form factors f B {p 2 ) and fr>{p 2 ) obtained with ([14]) 
and ( |T5| ) is shown in Fig. la and Fig. lb, respectively. Whereas the direct estimate of the 
form factors such as (P) is trustful at low and intermediate values of p 2 , but breaks down 
as p 2 —>■ ml and m 2 , respectively, there are no convincing arguments in favour of the pole 
model to be valid at p 2 — > 0. Hence we suggest to match the two descriptions in the region 
of intermediate momentum transfer near the upper end of the full curves plotted in Fig. 1. 
More precisely, in the following applications we use the direct calculation of f B (f B ) up to 
p 2 = 15 (1.0) GeV 2 and the single-pole approximation at higher values of p 2 . 

In order to improve these calculations further one has to take into account the 0(a s ) hard 
gluon exchanges in the correlation function (||]) explicitly. Parts of this not quite straightfor- 
ward task have already been completed [|14]]. 



3 Extraction of V uu from B — ► ttIu 

With the form factor f B {p 2 ) a t hand, we are now in the position to determine the CKM 



parameter \V u b\ from the recent CLEO measurement of B° 
charged lepton energy spectrum in the B° rest frame, 



7r Fig. 3a shows the 



dT(B° -> n-l+i/i) G 2 \V, 



ub 



dE, 



lQir 3 m B Jo 



dp 2 



2E l (m 2 B -ml + p 2 



-m B (p 2 + AEf 



IfUp 2 



(17) 
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with Pmax(Ei) = 2Ei( m B — Tn^/{rriB — 22^)). Integration of the spectrum over E\ yields the 
exclusive decay width 

1X5° -> n-l+Vi) = 8.1 \V ub \ 2 ps- 1 . (18) 

Experimentally, combining the preliminary CLEO result, BR(B° — > tt~1 + vi) = (1.63 ± 
0.46 ± 0.34) ■ 10~ 4 §, with the world average of the B° lifetime [0, t b o = 1.57 ± 0.05 ps, 
one obtains 

T(B° -> 7T-/+Z//) = (1.04 ± 0.37) ■ 10~ 4 ps" 1 , (19) 

where the errors have been added in quadrature. We take the CLEO result obtained by using 
the BSW model || for the efficiency estimate since the BSW form factors agree quite well 
with our predictions. 

Comparison of fllSP with (|T9| ) yields 

|Kfe| = 0.0036 ±0.0006 . (20) 

The theoretical uncertainty which we estimate conservatively to be less than 20% is not in- 
cluded in ©. 

dT/dE l [| V ub | 2 -10- 12 ] 

5 



3 
2 
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Fig. 3. The charged lepton spectra in (a) B° — ► ir l + u and (b) D° — > ir l + u . 

In the case of the Cabibbo suppressed D meson decay D° — ► tt~1 + vi we predict the charged 
lepton energy spectrum shown in Fig. 3b and the integrated semileptonic width 

T(D° -> 7r~/ + ^) = 0.156 |Kd| 2 ps _1 = (7.6 ± 0.2) • 10~ 3 ps" 1 , (21) 

where we have substituted |V^| = 0.221 ± 0.003 [0. 



This prediction should be compared with the experimental result 

T(D° -> 7r-e + z/ e ) = (9.4^|) • 10~ 3 ps~\ (22) 

derived from the branching ratio BR(D° — > it~e + v) = (3.9l^|) ' 10~ 3 and the lifetime r^o = 
0.415 ± 0.004 ps [OJ. The present theoretical uncertainty in (|21|) is of the order of the 
experimental error in (|22|). 

We see that the CKM-suppressed exclusive semileptonic widths are not yet measured pre- 
cisely enough to really challenge theory. In our approach, when better data become available, 
one may reduce the remaining uncertainties considerably and determine the CKM matrix el- 
ements with better accuracy. The idea is simple. Both the B and D meson sum rules are 
determined by one and the same set of nonperturbative parameters, including the coefficients 
of the nonasymptotic terms in the light-cone wave functions as exemplified in (Jfty. The present 
uncertainty in these parameters may be further reduced by constraining them through a pre- 
cise measurement of D — > nlu. After evolution from \i ~ m c to /i ~ m&, one should then 
obtain a more accurate prediction on B — > irlv and thus improve also the result on \V u b\. 

In conclusion we mention that the second form factor f~(p 2 ) of the B — > 7r transition 



matrix element (0) can be calculated by the same method. The work is in progress [0. This 
form factor may become measurable in future experiments in the decay B — > titv t . 



4 The nonfactorizable B —* J/ijjK amplitude 

To proceed from semileptonic to nonleptonic decays when QCD is switched off, one simply 
has to replace the lepton current by the relevant quark current in the weak Hamiltonian 
and multiply the latter by the corresponding CKM matrix element. Of course, the strong 
interaction effects then become much more complicated. They arise from (a) hard gluon 
exchange at short distances, (b) soft interactions of quarks and gluons including nonspectator 
effects, (c) the confinement mechanism leading to the formation of hadrons, and (d) final state 
interactions among the hadronic decay products. 

Theoretically, only the effect (a) can be systematically taken into account in the framework 
of improved QCD perturbation theory. The result is an effective weak Hamiltonian at the 
physical scale fi << mw of interest, given by a sum of local operators with renormalized 
Wilson coefficients |18[ . In the following we consider the decay B — > J/ipK as an important 
and instructive study case. The part of the effective Hamiltonian relevant for this decay may 
be written in the form 

H w = ^=V cb V c * s {(c 2 + |)0 2 + 2 Cl 2 } , (23) 

where 

O a (Ai) = (cTPc)(sT p b), OM = (cT^c)(sT p ^-b) (24) 
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with T p = 7 p (l — 75) and A a being the usual SU(3) colour matrices. The Wilson coefficients 
Cj(//) contain the QCD corrections from short distances below a characteristic distance scale 
of the order of the inverse 6-quark mass. The remaining task is then to calculate the hadronic 
matrix elements of the four-quark operators (|23|) incorporating the effects (b) to (d). This 
clearly requires nonperturbative methods whence progress has been slow and difficult. 

In the most radical approach, the matrix elements of Hw are factorized into products of 
hadronic matrix elements of the currents that compose H^. Strong interactions at scales 
lower than \x between quarks entering different currents as well as nonspectator effects are 
completely neglected in this approximation. Hence, the resemblance of nonleptonic decays to 
semileptonic decays goes quite far. Formally, the matrix element of the operator 2 vanishes 
because of colour conservation, and the amplitude for B(p + q) —>■ J/ip(p) K(q) (momenta in 
parenthesis) is given by the factorized matrix element of the operator 2 : 

(J/^K \H W \B) = ^=V cb V c * s ( c 2 (/i) + ^] (J/ijK I 2 (jjl) I B) (25) 



with 

{J/i)K I 2 {fi) I B) = (J /if; I cT p c I 0)(K \ sT p b | B) 

= 2Uf^ K m^-q) . (26) 

While the decay constant can be obtained from the leptonic width T(J/ip—* the 
B — > K form factor fg^ K (p 2 = mL) can be calculated in analogy to the B — > % form factor 
discussed in Section 2 . The result used below is fB^K( m %) = 0-55 ± 0.05 Q. Obviously, 
denotes the J/ip polarization vector. At this point one encounters a principal problem: 
since the matrix elements of quark currents in ( ^B]) are scale independent, the /^-dependence 
of (J/tpK I 02(/i) I B) which may cancel the /^-dependence of the Wilson coefficients in 



(|25|) and give a physically sensible result, is lost (see also the discussion in [0)- Hence the 
approximation fl2"6j ) can at best be valid at a particular value of fi which could be called the 
factorization scale \ip. Usually, /xp is expected to be of the order of the heavy quark mass. 
Using the next-to- leading coefficients c\^{^) in the HV scheme for = 225 MeV from 



[|19f| and taking /t = m& ~ 5 GeV, one obtains a 2 = 0.155 and from that the branching ratio 

BR(B -> J/'ipK) = 0.025%. (27) 



This estimate is significantly smaller than the measured branching ratios jib], |2"U|: 



BR(B~ -> J/ipK-) = (0.102 ± 0.014)% , (28) 

BR(B° -> J/i)K°) = (0.075 ± 0.021)% . (29) 

This failure, together with the scale problem pointed out above, indicates that factorization of 
matrix elements has to be accompanied by a reinterpretation of the Wilson coefficients. For 
class II decays such as B — > J/tpK, the short- distance /i-dependent coefficient c 2 (fj) +ci(/x)/3 
has to be substituted by an effective coefficient a 2 . Phenomenologically 0, the latter is treated 
as a free parameter to be determined from experiment. From fl2"5]), (PB"D, and (P%|), the most 
precise of the two measurements, one finds 

\a 2 \ = 0.32 ±0.02 (30) 

where the quoted error is purely experimental. The sign of a 2 , of course, remains undetermined. 
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More theoretically, one can argue that the above findings are a clear indication for the 
existence of significant nonfactorizable contributions which account for the difference of a 2 
from the value of c 2 + Ci/3, by roughly a factor of two at n — m b , and which soften the 
strong /i-dependence. A deeper analysis shows that the dominant effects should arise from the 
nonfactorizable matrix element of the operator 2 . Adopting the parameterization 

(J/^K | 2 (/i) | B) = 2f(v)Um^ ■ q) , (31) 

one derives from (|23|) and (^6]) the effective coefficient 



a 2 = c 2 (//) + — — + 2ci(/i)— — . (32) 

JB^K 

We consider the proof of this conjecture and the direct calculation of a 2 (and of analogous 
coefficients relevant for the other classes of two-body decays |§) as one of the most important 
tasks in heavy flavour physics. As a first step in this direction, we have undertaken a rough 



estimate of / using again QCD sum rules [21|. Following the general idea put forward in 



we work with the four-point correlation function 

< | T{j«(x)jt(y)6 2 (z)j*(0)} | > , (33) 

where jf£ 5 = uj^^s , = c^ u c, and jf = bi^u are the generating currents of the mesons 
involved. 

On the one side, one considers a dispersion relation for this function in terms of intermediate 
hadronic states in all three independent kinematical variables (the squared momenta of the 
currents). Similarly as in (§), the ground state contribution contains the matrix element of 
interest, that is (J/tpK | 2 (/i) | B). In addition, one has contributions from all kind of 
excited resonances and continuum states with a very complicated singularity structure. 

On the other side, in the deep-euclidean region one can apply a short-distance OPE to 
the correlation function (|33]). This leads to a representation by a series of vacuum expec- 
tation values (VEV) of local operators with calculable coefficients. The VEV are universal 
nonperturbative parameters characterizing the QCD vacuum. Best known examples are the 
gluon and quark condensates. We have performed this calculation including all operators up 



to dimension 6. More details of our analysis are discussed in [21]. Here we only mention two 
complications that are not present in the more familiar two- and three-point sum rules. One 
problem is the presence of a light continuum in the B channel below the 5-meson pole. This 
contribution can be associated with processes of the type i D*D s 7 — > J/tpK where an initial 
state carrying D*D S quantum numbers which is created by the combined action of the opera- 
tor product 2 ii from the vacuum (in (|33| ) ), rescatters into the final J/ipK state by strong 
interaction. As a reasonable solution we suggest to cancel this unwanted piece against a corre- 
sponding term in the OPE of (|33|) which develops a nonzero imaginary part at (p + q) 2 > Am 2 
in the channel with the corresponding quark content ccsq. The second problem concerns the 
subtraction of the contributions from higher resonances and continuum states in the remaining 
sum rule. Since the normal procedure (mentioned in section 2 ) does not work here, we use 
a rough model for the spectral function including besides the B, J/tp and K ground states, 
only the first excited resonances. We then perform a Borel transformation in the I?-meson 
channel and take moments in the charmonium channel. In the i^-meson channel, we keep the 
momentum variable q 2 spacelike. We then fit the spectral function to the sum rule varying 
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the Borel mass M, the moments and q 2 . From this fit we find / = — (0.045 -j- 0.075), where the 



implicit scale \i is identified with the central value M ~ yf^B ~ m b — 2-5 GeV. Substituting 
this estimate in ([32]), and evaluating the short- distance coefficients ci^/x) also at /i = M 



instead of the higher scale /i = mj used in fl27|) , we obtain 

a 2 = -0.29 + 0.38 - (0.19 4 0.31) = -(0.10 4- 0.22) , (34) 



where the three terms refer to the corresponding terms in (|32|). 

Several comments are in order. Firstly, the nonfactorizable matrix element (j31[) is con- 
siderably smaller than the factorizable one given in ( pE|) , numerically, |///B_»i<r("T.^)| — 
0.1. Nevertheless, it has a strong quantitative impact on a 2 because of the large coefficient 
|2ci/(c2 + ci/3) | ~ 20 -T- 30. Secondly, the factorizable term ci/3, nonleading in 1/N C , and the 
nonfactorizable term proportional to / are opposite in sign and hence tend to cancel. In fact, 
if |/| is taken at the upper end of the predicted range, the cancellation is almost complete and 
as a result increases the branching ratio considerably. This is exactly the scenario anticipated 
by the 1/A^-rule |23j and supported by similar estimates of nonfactorizable amplitudes for 



-D-decays |22[ and for B — > Dn |^4[ [25| . Thirdly, the theoretical value ( p4[) of |ct2| is still some- 
what low as compared to the phenomenological value (|3^). One should however note that the 
improvement due to the nonfactorizable contribution is remarkable as the factorizable terms 
alone would give | a 2 1 = 0.09 for // = M. Finally, our estimate yields a negative overall sign for 
a 2 in contradiction to strict factorization (first two terms in ([52])) and also to a global fit of 



the data |2(J, We stress, however, that in this fit the positive sign of a 2 is essentially deter- 
mined from the channels B~ — ► D°tt~ , D°p~ , D* tt~, and D*°p~ and is then assumed to hold 
universally for all channels. This assumption may not be correct. Certainly, the theoretical 
approach described here provides no justification to expect universal values and/or universal 
signs for the coefficients cij in different channels. From the relation (|32|) we see that univer- 
sality can at most be expected within certain classes of decay modes, separately. This point 
of view is also supported by the phenomenological analysis of nonfactorizable contributions in 



23 



As a last remark we emphasize that there is no simple relation between B and D decays 
in our approach since the OPE for the corresponding correlation functions significantly differ 
in the relevant diagrams and in the hierarchy of mass scales. We hope to clarify these issues 
further. 



5 Conclusion 

QCD sum rule techniques are very useful in calculating hadronic matrix elements for exclusive 
heavy meson decays. In particular, in combination with light-cone wave functions for pions and 
kaons they provide a powerful tool to derive heavy-to-light form factors and strong couplings 
of heavy to light mesons. In contrast to other methods including HQET, one can obtain 
predictions for the whole range in momentum transfer, and also for D mesons which are 
presumably too light for HQET to work. Moreover, finite mass effects are automatically 
included. 

In this paper, we have focussed on applications to the CKM-suppressed semileptonic decays 
B — > nlui and D — > nlvi, and to the factorization approximation of the nonleptonic decay 
B -> J/ipK. 
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Using the latter mode as a prototype example, we have discussed the shortcomings and 
principal problems of factorization of nonleptonic amplitudes. Again with the help of QCD sum 
rule methods, we have directly estimated the nonfactorizable contributions to the amplitude 
for B — > J/ipK. Our results indicate that nonfactorizable effects indeed play an essential role, 
putting the l/iV c -rule to work similarly as in D-decays. Moreover, on theoretical grounds the 
nonfactorizable amplitudes are expected to be channel-dependent. Hence, the pattern of two- 
body weak decays may be much richer than what is revealed by the current phenomenological 
analysis of the data. 

Acknowledgements. R.R. thanks Neville Harnew and Peter Schlein for the opportunity 
to participate in this fruitful and enjoyable workshop. 
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